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PREFACE 

nila  Memorandum  presents  the  results  of  a  theoretical  study  of  the 
behavior  of  certain  types  of  transmitting  antennas.  The  material  has 
possible  eventual  application  to  space  casssmleatloiis  and  to  problems  la 
eosnsmloatloBS  eaglaserlng  and  selenoe. 
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Tills  Memaraaftun  deals  with  the  deterolaatlon  of  the  prohablllty 
distribution  of  the  electric  field  resulting  from  an  arbltrazy  random  array 
of  sources  (or  scatterers)*  The  distribution  Is  surprisingly  simple.  Is 
easily  calculated  for  most  Interesting  eurray  distributions,  and  has  vide 
generality  of  application.  Specifically,  ve  find  the  antenna  pattern 
distribution  of  a  synthetic  aperture  antenna  fomed  by  a  moving  space 
vehicle  emitting  pulses  randomly  In  time.  However,  our  results  apply  not 
only  to  synthetic  aperture  antennas  of  arbitrary  distribution  but  also  to 
randomly  deleted  antennas  and  to  chaff,  meteor  trail,  and  electron  cloud 
diagnostics  as  veil. 

TSie  problem  Is  restricted  to  the  study  of  the  far  field  from  H  sources, 
the  positions  of  vhlch  are  Indepenlent  Identically  distributed  as  7(r). 
Nsrkov's  method  Is  then  used  to  analyze  vhat  Is  essentially  a  two-dimensional 
random  valk  Induced  by  a  three-dimensional  distribution.  It  Is  shown  that 
If  the  Fourier  transform  t(lc)  of  the  distribution  function  F(r)  can  be 
performed  In  closed  form,  then  the  limiting  form  of  the  probability  density 
of  the  resultant  electric  field  vector  Is  lamedlately  obvious  for  every 
frequency  and  direction  of  propagation.  Finally,  the  probability  density  of 
the  resultant  power  or  envelope  Is  determined  In  closed  form,  and  the 
correlation  betveen  the  resultant  field  at  different  angles  and  fTeqiueneles 
Is  exhibited. 
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I«  JJHIPBUOTIPW 

This  Maioraad.um  Is  eoneeraed  vlth  the  probability  distribution  of  the 
far  field  resulting  fron  an  array  of  H  sources,  the  positions  of  vhleh  are 
Independent,  Identically  distributed  random  variables,  nie  Importance  or 
utility  of  this  sliqply  stated  problem  Is  essentially  derived  trm  the 
extresM  simplicity  of  form  of  the  answer. 

Naffett^^^  briefly  discusses  a  specific  random  antenna— a  linear  array 
with  H  sources  unlfoimly  distributed.  However,  the  analysis  Is  Incomplete 
In  that  he  derives  the  distribution  of  only  one  of  the  two  quadrature 
coeponents  of  the  field.  Hence,  nothing  can  be  said  of  the  distribution  of 
the  envelope  or  power.  In  this  Memorandum  an  ei^esslon  will  be  derived 
for  the  envelope  distribution  of  the  resultant  field  tram  arbitrary  spatial 
distribution  functions. 

An  Interesting  and  very  general  study  of  the  problem  has  been  made  by 

(2) 

Kelly  and  Lemer'  '  In  which  they  e3q>lolt  the  mathematical  analogy  to  the 
shot  effect  of  the  radar  echo  from  a  random  collection  of  scatterers.  How¬ 
ever,  their  specializations  are  made  to  oases  In  vhleh  the  expected  axxtexmn 
pattern  Is  essentially  zero. 

dhe  following  exaaqples  will  lUxistrate  some  of  the  eurray  distributions 
which  may  be  of  Interest.  A  gas  diffusing  from  a  point  source  might  have 
a  three-dimensional  gausslem  distribution.  A  synthetic  aperture  antenna 
formed  by  emitting  pulses  at  Poisson  Increments  In  time  from  a  linearly 
moving  space  vehicle  ccm  be  considered,  under  certain  assumptions  on  the 
data  processing  scheme,  to  have  sources  having  a  vmiform  probability  density 
over  a  line  segment.  A  short  section  of  meteor  trail  may  have  Its  scatterers 
obeying  a  \mlform  distribution  over  a  cylinder  cut  out  by  the  meteor.  Also, 
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a  phased  eorray  whose  conopoaents  are  razilaimly  deleted  or  deunaged  is  a  nmdon 
array  with  a  distribution  corresponding  to  the  related  conrplete  array.  !Ihe 

(3) 

esqweted  changes  in  gain  of  randomly  deleted  arrays  are  discussed  by  Ogg.' 

The  far  field  assinoptloni  made  in  order  to  give  the  axunrer  its  sim¬ 
plicity  and  utility,  restricts  some  of  the  above  examples  to  relatively 
unlnterestlzig  cases.  In  order  for  the  far  field  assunptlon  to  hold,  the 
lateral  extent  of  the  target  must  be  less  than  the  diameter  of  a  Fresnel 
zone.  For  exanple,  the  lateral  extent  of  a  gas  cloud,  meteor  trail  segment, 
or  synthetic  aperture  antenna— for  distances  on  the  order  of  1000  km  at 
wavelengths  on  the  order  of  1  m— must  be  less  than  1  km.  Then  again,  the 
qualitatively  interesting  cases  occur  when  the  extent  of  the  array  is  a 
small  number  of  wavelengths,  because  the  field  components  will  tend  to 
add  coherently  rather  than  incoherently  as  the  extent  of  the  eurray  decresuies 
Figure  1  Illustrates  a  typical  sample  random  aurray  where  k  Is  the 
direction  of  propagation  vector  with  magnitude 


and  points  in  the  direction  of  propagation.  Under  the  far  field  assunptlon, 
for  identically  polarized  sources,  the  resultant  electric  field  I  at  a  point 
along  the  direction  k  can  be  written  as  the  complex  nuriber 
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Fig.  I —  Sketch  of  array  ond  electric  field  component 
orieing  from  source  ot  position  £ 


For  radar  echoes  k  becoaes  2k  because  of  the  doubled  jatb*  For  a 
distribution  g(r) 
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n,  TBE  ragp  Dia'miHffPMi  wwm  a  haibom  aruki 

nil*  section  vlll  provide  the  details  of  the  derivation  of  the 
asyiqptotle  prohablllty  distribution  of  the  resultant  electric  field.  For 
purposes  of  exposition  «e  assune  1  Isotropic  radiators  of  equal  asvUtude^ 
operating  In  phase.  We  nay  relax  these  assuoptlons  later. 

Hie  very  general  oethod  of  characteristic  functions  vlll  be  used  here, 
l.e.f  the  characteristic  function  of  the  sum  of  Independent  (not  necessarily 
Identically  distributed)  random  variables  (in  this  ease  vector-valued  random 
variables)  Is  the  product  of  the  individual  characteristic  functions.  Hie 
Inversion  foniula'  '  yields  the  desired  distribution  function  of  the  vector- 

(5) 

valued  random  variable.  Chandrasekhar^ ascribes  to  Markov  the  origin  of 
this  technique  vlth  respect  to  multidimensional  random  variables. 

Hie  essential  slsqpUflcatlons  occur  when  the  N  random  variables  are 
Identically  distributed  and  H  Is  larcpi  (lO  or  more)^  thus  allowing  an 
asynqptotle  expansion  of  the  characteristic  function^  vhlch  Is  then  easily 
Inverted. 

Our  problem  Is  essentially  tbat  of  solving  a  two-dimensional  random 
walk  (each  step  being  an  electric  field  vector)  Induced  by  the  three- 
dimensional  probability  distribution  of  the  positions  of  the  radiators.  We 
assume  the  use  of  Isotropically  radiating  sources  of  equal  Intensity  that 
are  operating  In  phase.  Hiese  sources  are  Independently  identically  distri¬ 
buted  In  space  according  to  the  probability  distribution  F(r)  where 
r  B  (x,y,s)^f  a  three-dlawnslonal  column  vector.* 

J<£)  -  Py  {x  <  X,  T  s  y,  Z  a  zj 


Hm  spidiol  "t”  designates  tzmnspose. 
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To  thli  dlatrlbutlen  fUBetlm  F(r)  tlunr*  eomspenda  the  eharaeterletlo 

fuxustlQD^^'^^  t(k)  ■  t(k„j  k  j  k  ).  If  the  radlatore  aleo  have  aone 
*  X  y  z 

randoameae  of  pbaae  6.  ve  nay  define  a  nev  k  ■  (k.,  k  .  k.,  l)^  and  a  nav 

—  X  y  8 

£  "  y*  >>  and  treat  thla  oaae  aleng  vlth  the  fenaar.  IB  either 
OBB#  we  hare  the  definition 

t(k)  .  X  {e  — }  -  (5) 


Ihe  contribution  to  the  electric  radiation  field  In  direction  k  due 

Ik^r  Ik^r 

to  a  aource  at  yoaltlon  r  la  e  —  — .  Ihe  eoaqplcx  nunber  e - vlU  be 

vrltten  aa  the  tve-dlaenalenal  column  vector  E  -  (coa  k^r,  aln  k^r)^. 

Hence  the  total  electric  field  vlU  be  the  aum  of  R  tvo-dlnenalonal  vectora^ 
each  of  unit  length  but  having  randcai  phaaea,  the  dlatrlbutlona  of  which 
depend  on  P(r)  and  k*  Our  flrat  taak  la  to  obtain  the  prohablUty  denalty 
p(l)  T»haro 

R 

s-iEii 

1^1 

Vote  that  |x|  ■  1  If  all  the  aoureea  add  coherently.  Eanee,  ve  call  X  the 
nenallaed  reaultant  vector. 

It  VlU  be  vlae  to  ceoalder  first  the  dlstrlhotlon  of  the  centered 
randen  variable 


R 


(7) 


6 


%  ■  ‘  {Si} 


iB  the  expected  value  of  the  electric  field  in  direction  k  due  to  the  1-th 
source  (averaged  over  the  spatial  distribution  F(r)).  Later,  a  change  in 
variable  will  yield  the  desired  distribution,  flhe  characteristic  function 
•(j)(w)  of  is  given  by 


CD 

«e» 

®  N 

^  J*1 


(8) 


vhere  the  independence  of  the  r^  has  been  used.  Making  use  of  the  identical 
distribution  of  the  r^  the  foUovlng  equation  is  true  for  any  Ji 


This  Integral  can  be  evaluated  only  for  very  special  fonss  of  P(r).  See 
Section  VI  for  such  a  case. 

In  order  to  proceed  vlth  the  analysis  for  general  F(r),  consider  the 
expansion  ofSq.  (9) 


CD 


(JO) 


7 


'baeauae  1 


{sW  • 


0.  Then 


f(,)(2)  -  exp  [-  *  0  (s-^^ 


vhleh  tcoda  In  the  Halt  for  large  H  to  the  characteristic  function 


r  2  SS-i 
*(w)  -  exp  “5“J 


vhere  Q  la  the  2x2  covariance  natolx  corresponding  to  a  single  radiator 


of  unit  a^Utudet 


Vhr(Ej^)  Cov(EpBj) 

Cov(BjEp)  Vhr(Kj) 


(incidentally,  the  asyngitotlc  derivations  of  results  similar  to  these  In 
Ref.  5  are  Incorrect.  )  The  characteristic  function  (Bq,  (U))  Is  easily 
Inverted^^^  to  give  the  density 

K2(i))  ■  1 2<»)  1(H)J 

A  change  of  varlahle  correspoodlng  to  definitions  found  In  Iqs.  (6)  and  (7) 


results  In 


*Tbe  asymptotic  e^nslons  are  Incorrect  in  Eqs.  (85),  (91),  and  (96), 
although  the  ansvers  exhibited  in  Eqs.  (87)  and  (93)  are  correct.  The  general 
answer,  £q.  (103),  is  Incorrect  because  Eq.  (lOO)  is  a  ncment  natrlx  rather 
than  a  covariance  matrix. 


8 


So>  for  large  the  probability  density  of  the  normalized  resultant 
electric  field  from  N  Independent  Identically  distributed  sources  Is 
biimrlate  noorml  vlth  sMsa  sad  oovarlaaee  matrix  ^  Q. 
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m.  miflATiOT  OF  yaua  aid  oo^gmiiuroB  vmax 


Ita  this  seetloa  the  results  are  pleaslag  in  their  slapllelty  sal 
Interesting  In  their  laberpretatlon*  We  vlsh  to  flai  jt  and  Q  In  tens 
of  F(r). 

ly  the  definition  found  la  Ig»  (8) 


■  |eos  S  j^eln 


-  i(je)  ■ 


where  t(k)  Is  Jixst  the  characteristic  function  of  p(r)  (see  Eq,.  (5))  «ad 
JH(k)  Is  the  tvo-dlnenslonal  vector  formed  from  the  real  and  imaginary  pairts 
of  f.  Ve  conclude,  vqpon  comgparlson  with  Eq.  (3)>  that  the  expected  antenna 
space  factor  tram  the  reuidom  eorray  with  distribution  F(r)  corresponds 
exactly  to  the  actual  space  factor  from  a  source  intensity  distribution 
g(r)  equal  to  the  probability  density  P^(r).  51ils  is  intuitively  obvious* 
The  covariance  matrix  Q  also  Is  easily  expressed  In  teims  of  j[(k). 
prom  Iq.  (12)  ve  have,  for  a  single  unit  radiator,  the  matrix 


Vi-w 


Vi-tt 


(15) 
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tff  oaleulAtlBg  for  exaiqple,  ve  flnA 


MOB 

OB 

■  "i  ^  ^  ^  ^  ^£)  ( 

■  5  * 

-  Ep^  -  I  +  ♦p(^)  -  tp^Ck)  (17) 

Froceedlag  In  the  saae  oemier,  ve  finally  obtain 


5  *  5  V^)  -  ♦/(») 


Q  - 


I  tj(a)  -  tjCk)  ti,(k) 


I  tjCa)  -  ti(k)  tpCk) 

5  -  5 


E  ■  CV-'* 


(Ifl) 


anl 

*<I)  -  »(e»  f  <i) 

vhere  N(^,  ^  Q)  la  the  normal  density  function  with  mean  ]i  and  covariance 
matrix  Q.  A  camnent  on  the  sliiQ>llclty  of  this  result  Is  In  order,  nie 
cbaracterlstlc  function  t(k)  just  the  three*<dlaenslonal  Fourier  trans* 
font  of  F(r)  and  Is  precisely  the  far  field  pattern  for  a  nomrandOB 


u 


lUunlBatlon  7(r)*  Eere^  ve  find  that  'both  ^  and  Q>  and  therefore  P(e)> 
are  aliqple  functions  of  t  evaluated  at  the  argument  k  which  "Just  happens" 
to  'be  the  direction  of  propagation  vector.  Thus^  If  the  Fourier  transform 
of  the  dlstrl'butlon  function  F(r)  ccm  be  performed  In  closed  form,  then  the 
probability  density  of  the  nonaallzed  resultant  electric  field  vector  Is 
lamwdlately  obvious  for  every  frequency  and  direction  of  propagation. 


UnyORM  LacAR  ARRAY 


We  are  eoacemed  vlth  the  space  factor  of  N  Isotropic  sources  operating 
In  phase  vhoee  positions  r^  are  Independent,  Identically  dlstrlhuted  randosi 
varlahles  unlfOansly  distributed  over  an  interval  of  length  L.  (See 


rig.  2.) 


Z  8(y)  6(*) 


0 

V 


a  X  a 


L 

2 


othervlse 


Actually  it  was  this  problem  (suggested  by  Irving  S»  Reed)  which  led 
to  the  Investigations  reported  In  this  Memorandum.  It  vould  be  of  Interest, 
for  exanple,  to  know  the  radiation  pattern  of  a  synthetic  aperture  antenna 
realized  by  transmitting  pulses  at  Poisson  Increments  In  time  from  a 
rapidly  orbiting  satellite,  deep  space  vehicle,  or  meteor  trail.  Ihe  data 
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•re  slnply  processed  as  If  the  pulse  posltloxu  vere  uoknoim.  Such  a  raniosi 
aoteima  will  have  the  prqperty  of  suppressing  the  grating  lobes  arising 
from  transmission  of  regularly  spaced  pulses* 

To  proceed,  ve  find  simply 


^  •  t 

t(k)  -  ♦(  ^  I  J  *  ”  "  ^ 


L 


sin  k  ^ 

X  2 

"T  IT 

^x? 


(20) 


where  k  ■>  —  cos  a  and  a  Is  the  angle  Included  between  the  direction  of 

X  X 

propagation  vector  and  the  line  of  sources.  We  are  motivated  now  to 
define 

k  L  . 

X  L  cos  or 

^  2«  X 


(21) 


and  call  It  the  "effective"  antenna  length  In  wavelengths. 


*<5)  -  »  (e>  5 


E 


sin  K 
rt  Y 

0 


Q 


“  2 
1  .  1  sin  2itv  Bin  jtv 


1  1  sin  2irv 

2  "  5  2iry 


(22) 


0 


14 


Observe  that  for  y  ■  1#  3»  ••• 


Figure  3  shevs  the  sketch  of  the  expected  value  of  E,  alotig  with  the  one 
standard  deviation  lines  of  the  real  and  imaginary  part  of  E. 


Remark  1:  Here  P(e)  depends  on  the  direction  of  propagation  k  through  y. 
Remark  2:  A  meteor  entering  the  atmosjAiere  emd  creating  an  Ionization  path 
will  result  In  a  Y  which  Increases  linearly  with  time,  nxus,  a 
receiver  will  observe  a  return  similar  to  Fig*  3  considered  as 


a  function  of  tine 
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aUDSSIAW  ARRAY  IH  THREE-DIMEHBIOKS 

It  Is  clear  that  echo  statistics  from  an  array  of  nonnally  distributed 
sources  (or  scatterers)  eire  of  more  than  academic  Interest.  Such  a  distri¬ 
bution  corlses  naturally  for  a  gas  diffusing  from  some  source  point,  dipoles 
subliming  from  an  orbiting  source,  and  In  "nice"  explosions  (relatively 
colUslonless  esqpanalon  of  a  gas  which  was  confined  and  In  equilibrium  at 
time  zero). 

Kre  we  consider  the  distribution 


("3|r|2/2a2)  (24) 

and  eully  calculate  by  Eq.  (^)  the  correspondliig  characteristic  function 

♦(k)  »  e:®  [-|k|®  (25) 

Since  |k|  =  ~  ,  we  cure  led  to  Introduce  ^  ,  the  cloud  "radius"  In 
wavelengths: 

Thus  i(k)  =  e“^ 

and  t(^)  = 


p(E)  -  H(ii,  Q) 
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(27) 


Hm  abovs  results  sre  plotted  in  Fig.  U, 


/i.  ,  Mtan  valus  of  E 


Fig.  4  —  Gaussian  array  (N>IO) 


Remark  1:  P(e)  Is  independent  of  direction  of  propagation  k. 

Remark  2:  A  g£is  of  N  particles  diffusing  from  an  initial  point  concentration 
produces  a  y  which  is  proportional  to  the  square  root  of  the 
elapsed  time.  Now  the  echo  statistics  can  be  read  from  the  grapli. 
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3t  For  Y  >  1»  ^  eovelope  r  of  th«  electric  field  E  la  eseentleUy 
E^rlelgli  distributed. 

/X  r  r  ^ 

*  HTTITT  L - TETin — J  ^0 


mttT' 

nils  la  turn,  for  y  »  1»  tecoaes  P(r)  ■  2Br  e* 


UHFORM  DISTRIHJTIOE  OVER  A  DISC 

nie  statistics  of  the  imlform  distribution  over  the  area  of  a  circle 
are  of  Interest,  for  exBaQ>le^  If  one  has  coherent  transmitters  scattered 
randomly  In  the  oceem,  desert,  or  on  the  surface  of  another  planet,  or 
If  one  has  the  analogous  problem  of  auilyzlng  echoes  from  clusters  of  points 
on  the  sea.  Consider  a  source,  the  position  of  which  Is  a  random  variable 
having  a  unlfom  distribution  over  a  disc  of  radius  R^. 


1  2  2  2 

;  z  .  0,  a  eJ 

«V 


F'(r) 


0 


otherwise 


♦(k) 


2  Ji(itY) 

jJy 


(28) 


(29) 


2  Rq  sin  0 

where  J2(z)  Is  the  Bessel  function  of  order  one;  y  - - - -  = 

"effective"  apertiure  diameter  (wavelengths);  and  0  Is  the  angle  between 
k  and  the  z  axis.  Equation  (29),  together  with  Eq.  (l8),  yields  the 
probability  density  of  the  field  as  a  function  of  y. 


IS 


CfEBER  DISTRIH3TI0MS 

Ab  bas  been  pointed  out,  the  only  difficulty  In  finding  P(e)  Is  In 
coBaiputlng  the  Fourier  transform  of  F(r) .  This  has  been  done  in  many 
Interesting  cases,  and  no  further  exsagiles  will  be  emaseratad. 
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V.  ENVELOPE  DIBTRIBUTION 

TOie  probability  distribution  of  the  electric  field  vector  is  useful 

when  conipletely  coherent  reception  (detection)  is  realized,  l.e.,  the 

lAiase  of  the  carrier  is  known  at  the  receiver.  In  general,  however,  the 

phase  of  the  return  Is  of  secondary  Importance  >  £uid  we  shall  be  most 

Interested  In  the  probability  distribution  of  the  envelope,  or  magnitude 

of  the  field.  Unfortunately,  averaging  out  the  phase  proves  to  be  a 

(8) 

difficult  unsolved  problem,'  *  where  satisfactory  esipresslons  exist  only 
for  special  cases.  In  this  section  we  attempt  to  determine  the  envelope 
density  p(r). 

Assume  for  simplicity  that  F(r)  is  symnetrlc  in  r;  that  is, 
dF(r)  ■  dF(-r).  (ihe  same  analysis  will  follow,  after  dlEigonalltation  of 
Q,  for  the  general  case.)  Ihen  t(k)  is  real  for  all  k  and 

(30) 


(31) 


(32) 


I  +  I  ♦(a)  -  ♦^(k) 


which  we  rewrite 


Ihen 


5  -  i 


^  0 


0  4- 


^  0 


We  bave  for  the  probability  density  of  the  nonallsed  electric  field 


vector 


P(I)  -  i  <)) 


175"®  *5 


sitCx^Xg) 


0^  - 

»1  >2 


Let 


>  r  cos  6 
-  r  sin  X 


Then  after  reduction 


P(3P#0)  “  ^  ®  ® 


-  A  cos  2d  +  B  cos  6 


vtaere 


C  - 


Hr 

(X1X2) 


175  exp 


N(Xi+Xg)r^  N  (k) 

^  H  h  2ip- 


HCXg  -  X3^)r^ 

^  “  — ITTT — 

^1  ^2 

N  r  i(k) 

X^  -  I  +  I  ♦(2k)  -  1r^(k) 
Xg  -  i 
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< 


I 


Mad 

OB 

♦(k)  ■  JjJ COB  k*r  dF(r) 

hOD 

In  order  to  proceed  In  the  determination  of  p(r),  ve  list  the 
following  pertinent  Bessel  function  relationships: 

Bessel's  Integrals  for  the  Bessel  coefficients  (Ref.  9,  p.  20) 


T  1  •  -  Iz  sin  * 


dft 


n  k  0 


n  >  0 


(36) 


The  series  definitions  of  the  coefficients  (Ref.  9,  pp.  15,  77): 


vm/1  _\n^2a 


■^n'*^  "  ^  ml  (n%  m)I 
iM) 

-  (-1)“ 


<»  /I 

f  \  r  '5  ^ 

^n'*'  "  mf  (Bfa)l 
I»bO 


■  VO 


V‘)  ■  ("0°  v*'^ 


f  n  2  0 


f  n  k  0 


n  k  0 


n  k  0 


(37) 


f 
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nift  Fourier  series  expensioiis: 


iz  sin  I 


■I 


z  cos  C  )  ^in  f  _ 
e  ■  /_,  e  I, 


(^) 


Then  the  marginal  density  of  the  envelope  is  given 


J2it 

O 


de 


1 

^  Jo 


-  A  COB  20  +  B  cos  0 
®  Clo 


By  Eq.  (38) 


>2jt 


Rearranging 


ByEq.  (36) 


B  COB  0 


d0 


l^eiOO 


p(r)  -  C  ^  I^(-A)  I2^(B)  (-1)* 


Patting  z  ■  -A  in  £q.  (37) 


P(r)  -  C  ^  1J,A)  1^{B) 


(38) 
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Plmlly,  "by  Eq.  (37) 

00 

p(r)  -  C  Iq(A)  Iq(B)  +  2  C  ^  ijA)  Ig^Cs)  (39) 

n^l 


The  expression  given  here  for  the  density  of  the  envelope  r  Is  the 
principal  result  of  this  section.  Althou^^  p(r)  Is  expressed  as  an 
Infinite  siim,  much  Infonnatlon  can  be  obtained  from  closer  examination. 

Consider,  for  exainple,  the  case  where  =  Xg  ®  l.e.,  the  real 

and  imaginary  conponents  of  the  electric  field  E,  ere  Independent  random 
verlables  with  Identlced.  variances  X  and  with  meeuis  i|r(k)  and  0,  respec¬ 
tively.  (Observe  that  this  case  occurs  for  the  field  from  a  random  erray 
only  when  t  (k).)  In  such  a  case  the  parameter  A  In  Eq.  (39)  Is 

equal  to  zero.  Since 


1.(0) 


n  *  0 
n  yf  0 


all  terns  but  the  first  drop  out  of  Eq.  (39)  leaving 


p(r)  -  C  Iq(B) 


for  r  i  0 


(UO) 


nils  Is  Just  the  Rayleigh  distribution  with  perameters  \|r(k)  and  x/K, 
which  Is  discussed,  for  exBsqple,  by  Rice  (Ref.  10,  p.  106) . 

Consider  next  the  case  where  the  meem  value  of  E  Is  zero,  l.e., 
t(k)  >  0.  Here  we  have  E  obeying  an  elliptical  bivariate  normal  law 
with  mean  zero,  corresponding  to  an  eiqiected  node  in  the  antenna  pattern. 


2k 


HOv  the  paraneter  B  in  Iq.  (39)  !■  mto)  all  tena  of  the  raa  but  tha 
first  are  zero.  !Ibua 


p(r)  -  C  Iq(A) 


Hr 


exp 


for  r  i  0 

* 


(1^1) 


Finally,  for  cooparlson  purposes,  we  append  em  approxlmtion  that  Is 
^1 

good  when  ~  is  on  the  order  of  one.  We  propose  to  reduce  the  elliptical 
gausslan  bivariate  distribution  of  E  to  an  "equivalent"  circular  distri¬ 
bution.  ibis  is  then  Integrated  directly  to  give  the  Rayleigh  distribution. 
There  is  no  restriction  to  diagonal  Q. 

A  locus  of  constant  probability  of  E  (see  Eq.  (13))  is  the  ellipse 
(E  -  4) '  Q”^  (E  -  4)  -  1.  Ihe  area  of  this  ellipse  is  it  where  Xj^ 

axil  Xg  the  first  and  second  eigenvalues  of  Q;  that  is,  and  Xg  axre 
the  lengths  of  the  major  and  minor  axes  of  the  ellipse. 

A  circle  of  the  sane  area  would,  have  radius  r  (the  harmonic 

mean  of  the  lengths  of  the  axes  of  the  ellipse).  Finally,  recall  that 
Xj^Xg  *■  ||Q||  ■  11^  simply  replace  the  elliptical 


(**2) 
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by  tlM  elreulur 


Bence  Eq.  (13)  becomes 


When  the  phase  Is  Integrated  out,  the  probability  density  of  the  envelope 
r  Is 


(45) 


VI.  AN  mCT  CASE 


In  the  preceding  sections  ve  have  been  forced  to  consider  a  random 

walk  with  unit  steps  having  an  eurbltrary  distribution  of  angle  at  each 

(9) 

step.  This  differs  from  the  classlcel.  case,  treated  by  Kluyver'^'  and 
Rayleigh, where  all  angles  are  equlprobable.  However,  ve  do  get  the 
Isoperlodlc  case  when  our  source  distribution  F(r)  satisfies  certain 
conditions  (roughly,  a  uniform  distribution  over  a  volume  whose  dimension 
along  the  direction  of  propagation  Is  an  Integral  number  of  wavelengths). 
To  Illustrate  this,  consider  a  reuidom  linear  array  where  N  sources  are 
distributed  Independently,  uniformly  over  a  length  L  (see  Fig.  2).  Ihen 
whenever  k  *  L  Is  an  Integer,  l.e.,  ^ ~  *«  n.  It  Is  apparent  that 

the  distribution  of  the  angle  of  the  unit  electric  field  vector  from  each 
source  Is  uniform  over  2ir.  Hence  Kluyver's^^^  results  apply  (see  also 
Ref.  10,  p.  242  and  Ref.  U),  and  the  probability  density  of  the  envelope 
r  a  |E|  Is  exactly 

v(.r)  -  r  X  JQ(rx)  jJ  (x)  dx 

nils  has  been  shown  to  approach  a  Rayleigh  distribution  In  the  limit  for 
large  E. 


27 


vn.  CORHELATIOir  IN  RADIATIOW  AT  DPTEREllT 
Aistjia  M)  PREqUENCIES 

The  knowledge  of  the  radiation  field  from  a  random  array  at  one 
direction  and  frequency  will  allow  an  Improved  estimate  of  the  radiation 
field  at  another  direction  or  wother  frequency.  It  is  intuitively  clear 
that  specification  of  the  radiation  field  in  a  given  direction  puts 
constraints  on  the  positions  of  the  sources.  Ihls  information  Induces  a 
conditional  prohablllty  distribution  on  the  random  field  at  new  cuigles 
emd  frequencies. 

Applications  of  the  generalizations  in  this  section  core  numerous. 

An  exEuqple  of  Interest  is  the  correlation  between  radar  returns  from  a 
cloud  of  fixed  random  scatterers  when  the  carrier  frequency  is  chcmged. 

An  "invisible"  cloud  at  one  frequency  will,  with  a  certain  probability, 
become  visible  at  another.  Let  us  also  consider  a  certain  scintillating 
target  model  conotposed  of  N  randomly  distributed  scatterers.  The  observation 
angle  will  change  as  this  target  moves  past,  and  we  might  wish  to  use  the 
past  information  in  the  optimal  way  in  order  to  predict  future  fading  and 
provide  for  a  good  tracking  capability.  In  the  former  exasple,  we  consider 
a  direction  of  propagation  vector  k  changing  in  magnitude  only;  in  the 
latter,  k  changes  only  in  cmgle.  We  shall  treat  both  cases  at  once. 

ISg#  •••  >  ^  l>e  m  direction  of  propagation  vectors  for  which 
we  should  like  to  specify  the  radiation  field  vectors  E(kj^),  E(^),  ...  , 

E(k  ).  Recall  that  |k.|  =  .  Then  for  large  N,  so  that  the  individual 

Xl 

E(k^)  are  each  approximately  nozaally  distributed  as  expressed  In  Iq.  (IB), 
it  is  obvious  that 

(^7) 
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itlMre  BOW 

a  -  aCJSe)#  •••  »  ad^O* 

and  Q  la  the  coaipoBlte  matrix 


«Vii> 


•nl  irture  ««eh  )  is  the  2x2  correlation  metrlx 

«(ii.ij)  -  *  {(5(ii)  -  uj)  (s(iSj)  -  Ej)*} 

and  Is  shown  eo^llcltly  In  Eqe  (3l)e  We  evaluate  the  upper  left-hand 
element  of  a  typical  conqponent  matrix 

Ou^ii'isj)  ■  *  {(Vii)  -  VisiOCVij)  -  Vij))} 

•  *  {Visi>  Visj>}  -  ♦R(ii)  VJsj> 


(W) 


(k9) 


W 

JjJ COS  (kj*r)  cos  (jfcj'l)  '*?(£) 


-  Vki)  Vkj) 


cos  +  i  cos  clF(r) 


(50) 
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"  2  5  "  '•'r^-1^  V-P 


Froeeedlag  slallarly  ve  obtain 


E(ii)  ■  (V-i^'  ♦i^iiO* 


!nxus,  the  general  Joint  distribution  of  the  E(k^)  has  been  demons tirated. 
From  Eqs.  (^7),  (^8),  and  (^l)  the  conditional  distributions  nay  straight- 
fhrvardly  be  obtained  (see,  for  exangple.  Ref.  6,  p.  31?)  • 
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VIII,  DISCUSSION 

The  problem  we  have  studied  in  this  paper  is  essentially  the  following. 
What  is  the  ito  field  from  an  array  of  N  (isotropic)  sources  (of  equal 
amplitude  and  operating  in  phase) ,  the  positions  of  which  are  independent. 
Identically  distributed  vector-valued  random  variables?  The  phrases  in 
parentheses  are  noncruclal  assumptions,  made  primarily  to  simplify  the 
exposition  and  are  easily  eliminated  for  greater  generality.  However, 
the  underlined  assumptions  are  essential  In  that  they  are  necessary  for 
the  simplicity  of  the  results. 

Isotropy  and  equal  amplitude  of  source  radiation  are  noncruclal 
assumptions  if  the  probability  distributions  on  these  parameters  are 
Independent  of  position.  As  was  mentioned  In  Section  II,  distributions  of 
the  phases  of  the  sources  as  a  function  of  position  may  be  admitted  by  a 
mere  reinterpretation  of  the  notation.  However,  without  the  far  field 
assumption,  p(E)  would  not  have  the  nice  explicit  dependence  on  ♦(k), 
although  it  would  still  be  gausslan. 

The  limiting  field  distribution  is  given  in  Eq.  (l8);  the  envelope 
distribution  in  Eq.  (39);  the  envelope  distribution  at  expected  nodes  of 
the  antenna  pattern  in  Eq.  (4l);  the  joint  distribution  of  the  field  in 
several  directions  in  Eqs.  (1^?)/  (*^8),  (51);  the  field  from  a  gaussian 
array  in  Fig.  4;  and  the  field  frco  a  linear  array  in  Pig.  3« 

In  summary  it  can  be  said  that  the  limiting  form  of  the  probability 
density  of  the  resultant  electric  field  vector  arising  from  an  array  of  N 
sources,  the  positions  of  which  are  independent,  identically  distributed 
random  variables,  is  bivariate  normal.  If  the  Fourier  transform  of  the 
distribution  function  F(r)  can  be  performed  in  closed  form,  then  the 
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probability  density  of  the  resultant  electric  field  vector  is  Immediately 
obvious  (Eq.  (l8))  for  every  frequency  and  direction  of  propagation. 

It  is  suggested  that  examination  of  antenna  statistics  in  the  near 
field  is  of  Interest  and  that  the  methods  of  Chandrasekhar  (Ref. 

Chap.  IV)  and  Kolmogorov  (Ref.  12,  p.  I7I)  mltpxt  be  of  use. 
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